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Abstract 

A directed animal is a percolation cluster in the directed site percolation model. The aim of 
this paper is to exhibit a strong relation between the problem of computing the generating 
function GF of directed animals on the square lattice, counted according to the area and 
the perimeter, and the problem of solving a system of quadratic equations involving unknown 
matrices. We present some solid evidence that some infinite explicit matrices, the fixed points 
of a rewriting like system are the natural solutions to this system of equations: some strong 
evidence is given that the problem of finding GF reduces to the problem of finding an eigenvector 
to an explicit infinite matrix. Similar properties are shown for other combinatorial questions 
concerning directed animals, and for different lattices. 

The author is partially supported by the ANR-08-BLAN-0190-04 A3. 

1 Introduction 

We are mainly interested in the study of the area and perimeter generating function GF Sq of 
directed animals on the square lattice Sq, but other lattices and questions will also be addressed. 
The computation of GF Sq is a central question in enumeration problems for directed animals on 
two dimensional lattices, since it is deeply related to the study of directed percolation on the square 
lattice. In this paper, even if we do not find an explicit formula for GF Sq , we show that to compute 
GF Sq it suffices to solve a quadratic system of equations involving 4 unknown finite matrices. We 
are unable to find a solution, but we provide some infinite size matrices which appear as the natural 
solution to this system of equations. They appear to be a fixed point of a rewriting system, the 
rewriting rules involving the tensorial product of matrices. We give strong evidence that finding a 
right and a left eigenvector to these matrices should lead to GF Sq . We hope that this gives some 
insight on the algebraic structure of this problem, and that this will allow some readers to compute 
GF Sq . 

In Section [6j we show that numerous similar problems can be treated similarly. 

The set of oriented graphs with no cycle and no multiple edges which have a finite or countable 
number of vertices and bounded degree is denoted Q. For any graph G = (V, E) in Q, V is the set 
of vertices and E C V 2 the set of oriented edges. The orientation of the edges leads to the notion 
of a descendant: for (xi,X2) E E, X2 is said to be a child of x\ and the set of children of x% is 
denoted Ch(xi). A directed path d in G is a sequence of vertices (xi, . . . ,Xk) such that for any 
/ > 2, xi E Ch(x;_i). The vertex x\ (resp. Xk) is called the origin (resp. the target) of d. 

Definition 1 Let G = (V, E) be in Q, and S be a subset of V . 

• A directed animal (DA) A with source S is a subset ofV containing S, such that for every a G A 



Figure 1: DA on the square lattice and on the cylinder Sq(3). The black dots are the cells and the white 
dots, the perimeter sites. The cylinder is defined as the square lattice with periodic conditions, meaning 
that the extreme vertical lines on the second picture are identified. The outlined cells constitute the 
source. 

there exists a directed path having target a and its origin in S entirely contained in A. The cardi- 
nality j^A of A is called the area of A. 

• A perimeter site c of a DA A with source S is an element of V \ A such that {c} U A is still a 
DA with source S. The set of perimeter sites of A is denoted P{A). 

We denote by Ag the set of finite DA on G with source S. The generating function (GF) GFg 
counts the DA with source S according to the area and perimeter: 

GPf(x,y):= £ x#V P(A) - 

Hence, the area generating function is GF|(i, 1). 

The search for a formula for GFg(x,y) may be seen as the combinatorial contribution to the 
study of directed percolation per site models. Indeed, on a probability space (Q,A, P) consider 
a random colouring of the vertices of V by the colours and 1. Formally, this is given by a 
family of i.i.d. Bernoulli random variables (B v (p),v £ V) indexed by the vertex set (we then 
have F(B v (p) = 1) = 1 — ¥(B v (p) = 0) = p). The directed percolation cluster with source 
v G V is the maximum DA with source S = {v} included in the set of 1-coloured vertices, that 
is {it £ y : B u {p) = 1} (the empty case, possible here, arises with probability 1 — p): denote 
it A?(p). Since for any DA A with source v, ¥(A v (p) = A) = p^ A {l — p)& p ( A \ the percolation 
cluster is finite with probability 1 if (1 — p) + X^ag^ g P^ A (^ —p)^ p ^ AS> = 1, which is equivalent to 

GF?i(p, 1 — p) = p. Hence a computation of GF^j(x,y) would probably allows one to compute 
the directed percolation threshold, and/or the associated critical exponent. 

Denote by Sq = (Vs q , -Esq) the directed square lattice where V$ q = Z 2 and 

^Sq = {{(xi,yi), ix 2 ,V2)) G Vs q , such that (x 2 ,y 2 ) - (xi,yi) G {(0,1), (1,0)}} . 

Surveys on the study of DA on two dimensional lattices exist: Bousquet-Melou [1] and Le Borgne 
& Marckert [15]. When S is reduced to a singleton, the area GF GF^ q (x,l) is well known, and 
numerous different approaches are possible to compute it: the gas approach (Dhar [91 HP] , but 
also [1] , [15] , Albenque pQ ) , heap of pieces approach (Viennot [18] ) , combinatorial decomposition 
(Corteel & al [7j, Betrema & Penaud [3]). On the other side, almost nothing is known about 
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GF S s q (x,y) (except for Bacher [2] who computed J2p(A)x^ a \ = §u ( x >^) 011 the S( l uare 

lattice with or without periodic conditions, proving conjectures by Conway E] and Le Borgne [14] ) . 
GF 5 s q (x,y) is not believed to be L>-finite. 

The aim of this paper is to use some algebra to search for a formula for GF^ q (x, y). We will use 
the idea of Nadal & al. [16] and Hakim and Nadal also used extensively in [1]. First, the work 
is done on a so-called cylinder Sq(n), a vertical strip of Sq with periodic conditions (see Figure [I]). 
Second, the corresponding for Sq is obtained by taking a formal limit since small DA on Sq(n) and 
Sq are the same. We will proceed similarly here, starting with GF Sq ^ n \x,y). 

In order to highlight the different considerations leading to the introduction of infinite matrices, 
we have decided to simultaneously treat the study of GF Sq ( n \x,y) and a case where infinite 
matrices can be avoided: the computation of GF Sq '"' (1, 1). This leads to a new derivation of 
GF Sq (x,l). 

1.1 Two gases 

Let G = (V, E) G Q be an oriented graph. Following the ideas developed in [15], we define 
two processes = (X v ,v G V) and Yq = (Y v ,v G V) indexed by the vertex set, and taking 
their values in {0, 1}. For this latter reason, the processes are called "gases", the value 1 (resp. 0) 
representing the presence (resp. absence) of a particle. 

Both processes and Yq are defined on a probability space (f2, «4,P), on which are de- 
fined some families of i.i.d. random variables (£ v ,v G V) indexed by the vertex set, where 
3 V = {B\{p) , B^iq)) is a pair of independent Bernoulli random variables whose parameters are 
p and q. 

Gas of type 1 : For any v £ V, set 

x v =bhp) n u-*c). a) 

ceCh(-u) 

That is, if X c = for all c G Ch(u), then X v = 1 with probability p; otherwise X v = 0. Since two 
neighbouring sites can not be simultaneously occupied, this model is called a hard particle model 
in the physics literature. 

Gas of type 2 : For any v £ V, set 

Y v = B\(p) min{Y c : c G Qh{v)} + (1 - B^B^q). (2) 

Here Y v is equal to min{Y^ : c G Ch(w)} with probability p, and to B^q). with probability 1 — p. 

Lemma 2 Let G G Q . If p G (0, 1) is small enough, both processes and Yq are almost surely 
well defined. 

Proof. We use the argument in [15] (the argument being already present in the PhD thesis 
of Le Borgne p3]). For both gases, when B\{p) = 1, the set of values {X c ,c G Ch(w)} (resp. 
{y c ,c G Ch(v)}) is needed to compute X v (resp. Y v ), but they are not needed when B\{p) = 0, in 
which case X v = and Y v = B^q). The fact that these "recursive definitions" ([!]) and ([2]) indeed 
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define some objects is not clear, but the values of both X v and Y v are certainly well defined if 
A v (p) is finite, since in this case, the recursive computation of X v (and Y v ) using the values of the 
children ends since the value of X and Y on perimeter sites of A v (p) - sites where B\ (p) is zero - 
is well defined. Hence, if the family of DA (A v (p), v E V) is a family of finite DA, both processes 
are defined. Now consider the standard problem of the directed percolation threshold. Let 

p crit = sup{p : P(W, \A v (p)\ < +00) = 1}. 

Since we assume that the maximum degree if the graph is finite, p cr it £ (0, 1]. For all < p < p cr u, 
X(j and Yg are then a.s. defined. □ 

A subset S of V is said to be free if for any si,S2 G S with s\ ^ S2, there does not exist 
any directed path with origin s\ and target S2 in G. The following Proposition says that the 
computation of the finite dimensional distribution of the gas of type 1 (resp. type 2) is equivalent 
to the computation of the DA GF according to the area (resp. area and perimeter), for general 
source. A DA is said to have over-source S if is a DA with source S' , with S' C S. If S' 7^ S, the 
set S \ S' is taken to be a subset of the perimeter. We denote by GFg the generating function of 
DA with over-source S. 

Proposition 3 For any directed graph G = (V, E) in Q and any free subset S of V , 

F{X v = l,veS) = (-l)# s GFf(- P ,l), 
P{Y v = l,veS) = GPf(p,(l- P )g), 

where the first equality holds for \p\ smaller than the radius of convergence o/GF<?(— p, 1) and the 
second one holds if < p < p cr it- 

Proof. The first assertion is proved in [15] (Theorem 2.7 for a single source, Proposition 2.16 for 
any source) on a general graph, and was already used in [3] on lattices for a single source. Dhar 
[9] , who made the connection on lattices between GF of DA and the problem of finding the density 
of a hard particle system on an associated graph, did not use the construction of the process X, 
but different considerations of the same process. The second assertion is well-known, and it is also 
proved in [T3] (Theorem 4.3) and valid on any graph of Q. The reason is simple: Y v = 1 if and only 
if B^iq) = 1 for all perimeter sites u of A v (p) (where v is considered as a perimeter site of A v (p) 
in the case where A v (p) is empty). □ 

Note 4 Although GF(x, 1) is a projection of GF(x,y), and GFg(x,y) can be computed easily 
thanks to (GFg,(x, y), S' C S), the gas of type 2 is not an extension of the gas of type 1. To 
compute GF G (p, 1) using the gas of type 2, q needs to be 1/(1 —p), which is larger than 1; this 
is not possible for probabilistic considerations. Nevertheless, given the polynomial form of T Y 
(see Formulas ^ and T Y still has a meaning when q = 1/(1 — p); for this value, it is 

no longer a positive kernel, but X^be-E n ^efb = 1 f or anv a ^ ^n- ^ non-negative solution p to 
[i = fiT Y still exists, as can be checked by rewriting the following system of equations: let, for any 
Cc{0,...,n-1}, 

W c = p({x £ E n :ieC =>Xi = 1}). 

p solves the system [i = /xT Y if and only if W := (Wc, C C {0, . . . , n — 1}) is solution to Wc = 
EdccCC 1 _ P)g)°^ D P D Wch(D), which when (1 - p)q = 1 is a rewriting ofW c = J2dcC P° W Ch(D)- 
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Clearly, W satisfies the same system of equations as GF^ := (GF^(p, 1),C C {0, ...,n— 1}). 
Since GF^< exists and is non negative, fi = /uT Y admits some solutions. The initial conditions 
W = GFg(p,l) = 1 allows one to identify the two set of series W and GFg. ,4s in Dhar JPJ/ 
or Bousquet-Melou J^]/, Wc has a product form on E n , meaning that (W a = a n Y\2=i Qa lt a i+1 mod n 
for some numbers Qo t o,Qo,i,Qi,o,Ql,l)- The behaviour of P(Y V = l,v G S) becomes singular at 
q = 1/(1 — p): some drastic simplifications of the involved algebra appear, but only for that value 
of q. This leads to a Markovian type structure of fj,. 

Some extensions of the gas of type 1 can also be introduced, for example X x = B p riceCh(a;)(^ ~~ 
X c ) + (1 — B p )Bq. These extensions are either of the same type as this gas, meaning that an easy- 
to-prove Markovian behaviour occurs, or will present the same kind of difficulties as for the gas of 
type 2 (as illustrated by the different cases discussed in Section^. We think that any gas allowing 
us to compute GF Sq '"'(x,|/) will be (at best) as difficult to describe as the gas of type 2. 

1.2 The square lattice : the cylinder approach 

We study some properties of the processes of type 1 and 2 defined on Sq(n), the square lattice 
with periodic conditions shown in Figure [TJ We first need to label the sites of this lattice and 
its rows, in order to describe its Markov chain structure, row by row: the Zth row is Row w = 
{(x mod n,l — x mod n),x G Z} and let x®(i) = (i mod n,l — i mod n) be the ith element of 
this row. The two children of x^(i) are x^ l+1 \i + 1) and x^ l+1 \i). Note that Row*-' +1 ^ is above 
Row^ in Figure [l} 

It is useful to provide a row decomposition of the processes Xs q ( n ) and Ys q ( n ). For brievity, let 
Z G {Xs q ( n ), Yg q ( n )} be one of these processes, and let Z^ l '(i) be the value of Z at x^(i), and 

Z (Z) = (z®(i),i = 0,...,n-l) 

the value of Z on Row^. Clearly Z^ depends on Z^ +1 ) and the values of the Bernoulli random 
variables v G Row^) only. Then, the sequence 2fl\ when I goes from +oo to — oo, is a Markov 
chain. The transition kernel can be expressed using ([I]) and Q : denote by E n the state space 
{0, l} n , by a n := (ao, • • • , a n -i), and b n = (6o, • • • , b ri -\) some generic elements of E n , and by ® 
the addition in Z/nZ. We have 

T z (b,a) := P(Z« = a | 7fl +v > =h) = \{ ij^, (3) 

and where T^ y z and T^ y z are explicit: 

Tf^ z = F(B 1 (p)(l-x)(l-y) = z) (4) 

= l(:r,y)^(0,0) 1-2=0 + l{x,y)=0 (P^z=l + (1 ~ P)^z=o) , (5) 

T? yz = F(B 1 (p)xy+(l-B 1 (p))B 2 (q)=z) (6) 
= lmin(x,y)=0 ((P + (1 ~ P)( l ~ Q))lz=0 + (1 ~ p)qU=l) 

+lmm(x,y)=l ((P + (1 " P)<l)lz=l + (1 " P)(l " ?)l*=o) • (7) 

If, for a given I G Z, Z(' +1 ) has distribution u, then the distribution /i of Z^ satisfies 

M (a)= £ i/(b)T£(b,a). (8) 



5 



For clarity, we use an operator-type notation and write \i = z^T z ; accordingly, /x = v (T Z ) K 
designates the law of Z(~ K ) if has distribution v. 
Using the row decomposition of Sq(n), we have 

Lemma 5 l)Letn>\ be fixed and Z E {Xs q ( n ), Ys q ( n )} . For p, q € (0, 1), the process (Z(~i\ j G 
Z) (indexed by decreasing j's) is a Markov chain with finite state spaces E n under its stationary 
distribution. Since this Markov chain is irreducible and aperiodic, the distribution of for any 
j is the only probability measure solution /2 Z of 

^ = ^l- (9) 

2) The distribution // z is the limit distribution of any Markov chain with transition T z on Sq(n), 
i.e. for any distribution v on E n , we have i/(T z ) K — > /U Z . 

The matrix T^, indexed by the elements of E n is called a (row) transfer matrix in statistical 
physics literature. By the Perron-Frobenius theorem, the equation ^ has a unique solution up 
to a multiplicative constant, which can be turned into a probability distribution. Finding a /x z 
solution of Q is a problem of linear algebra, and a solution can be found by using a computer for 
small values of n. The fact that T z has the product representation given in ^ plays a secondary 
role in that respect. In the following, it will play a role of primary importance. 

Proof. The only unclear assertion in (1) is that the row process Z^ - ^ follows under its stationary 
distribution. This comes from the infinite construction that we have, under which Z^ and Z^ +1 ^ 
have the same distribution. Assertion (2) is a well known property for aperiodic irreducible Markov 
chains on a finite state space. □ 

Note 6 In the case of the square (or triangular) lattice the gas of type 1 is a Markov chain on the 
horizontal lines of the lattice J75] /. In ^j, it is observed that it is a Markovian field on a zigzag 
formed with two consecutive lines ofSq(n); this Markovian field converges when the row size goes 
to +oo to a Markov chain on the line according to JI^. It turns out that the gas of type 2 is not 
Markovian on a line of the cylinder, on a line of the lattice, or on the zig-zag %4%. In other words, 
there does not exist any 2x2 matrix A such that 

n 

Hj(x!,...,x n ) = Cn]jA XitX ^ +1 . (10) 

i=l 

To show this, one can solve |9p with a computer for n = 3, 4, . . . , to check that no factorisation of 



the invariant measure compatible with (10) is possible (a similar work can be done directly on the 
rows of the entire lattice). It can also be checked that no solution corresponding to a Markov chain 
with longer memory exists, at least for small memories. When a computer is used to compute the 
invariant distribution on a cylinder of small size, no "regularity" of the measure has been observed: 
it turns out that on the cylinder of size 12, if fi^ixi, ■ ■ ■ , x n ) = /i^(yi, • • • , y n ), then x and y are 
equal up to a rotation and/ or symmetry. This implies that what happens is drastically more complex 
that what appears for the gas of type 1 where has the form (10). 



We now look at new considerations. 
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2 A new paradigm 

PropositionKHsays that /^f(l,*, ...,*) := Yj Xi e{0,i},i=2,..., n (!) x 2, ■ ■ - ,x n ), i.e. the probability 



that Y" = 1 at some position v on Sq(n) is, up to change of variables, equal to GF s ^ n \x,y). A 
nice description of would thus be helpful. For this, two new ideas arise: 

■ The first one is to search for a representation of the same type as (10), but for matrices A x ,y 



(taking the trace afterwards) instead of "real numbers A xy " . This is not just a way to add some 
degrees of freedom: Proposition [7] below says that all measures invariant under rotation on E n have 
such a representation, and moreover have a representation of the form Trace(f3^ 1 Q Xi ) for some 
matrices Q° and Q 1 , 

■ The second idea comes from the product form ^ of the transition kernel for Z G {X, Y}. 
It suggests that a "local equation" linking the matrix Q x "weighting the probability to observe x 
somewhere on Row(Z)" and Q y , Q y , the matrices weighting their two children in Row(Z + 1) could 
suffice. But two neighbours in Row(Z) share a neighbour in Row(Z + 1), and this must be taken 
into account. The second idea is to share the neighbour, its suffices to "split" the matrices, and 
search for a solution of the form Q x = V X H X for some matrices V x and H x having right splitting 
properties. This is what is done and proved to be possible, up to taking infinite matrices Q x . 

2.1 Transfer matrix and measure representations 

We recall the definition and some properties of the Kronecker product (or tensor product) 
between matrices that we will use in the paper. If A is an m x k matrix and B is a p x q, then the 
Kronecker product A B is the mp x kq matrix 



A®B 



AiiB ••• A hk B 



Am } iB ■ ■ ■ A m ] t B 
The Kronecker product is associative. For matrices A, B,C, D, 



{A®B){C®D) = {AC)®{BD), (11) 
Trace(A ® B) = Trace(A) Trace (S) (12) 



where it is assumed in ( 11 ) that A, B, C, D have sizes such that AC and BD are well defined, and 
in (12) that the matrices are square matrices. The Kronecker product extends to infinite matrices 
A ® B, but it is "interesting" only for finite B. 

For n > 1, let M n be the set of probability measures on E n , invariant under rotation: if 
H G M n , then for any x = (x , . . . ,x n -i) G E n , fj,(x) = fi(xi, . . .,x n -i,x ). 

Proposition 7 For any fi G M n , there exists two square matrices Q° and Q 1 of finite size, and 
four rectangular finite matrices V , V , H , H 1 such that: 
(1) for any x G E n , 

/n-1 \ 



/x(x) = Trace JJ Q Xi 



vi=0 



(2) Q 1 = V 1 !! 1 , Q° = V°H°, V 1 ^ = [o] and V^H 1 = [o] , where [0] stands for the null matrix 
with the appropriate size. 
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Proof. First, assume that there exists some matrices (Q x ,x £ {0,1}) of size m x m, such that 
//(x) = Trace(Ui=o Q Xi )- Take V 1 = Q 1 ® [l 0] , H 1 = Id(m) , U° = Q° ® [0 l] and 



H° = Id(m) 



where Id(m) is the identity matrix of size m x m. Using Trace(J} a« (8) bi) = 
Trace(n ctj) Trace(J| hi), we see that for any x E £" n , 

Trace (jj Q x ^j = Trace \JJ Q x ^j , 

and clearly Q 1 = V l H l , Q° = V°H°, V Q H l = V l H° = [0] . 

Hence, only (1) remains to be proved. For this, denote by M' n the set of probability distributions 
having the representation /i(x) = Trace(f3™ = g 1 Q Xi ) for some finite matrices Q° and Q 1 . We show 
that M.' n = M. n . For this, observe that M! n is closed under finite mixture: if \i = a^ a + (1 — 
and /Xj(x) = Trace(n"r o Q Xl ) for j 6 {a, b}, then taking 



a^Q* 
P 1/n Q x . 



XG{0,1}, 



we obtain fi(x) = Trace(n™_T Q Xi ). This shows that Ai' n is closed under finite mixture. 

Since Trace(^4-B) = Trace(-EM), any probability distribution in A4' n is invariant under rotation. 
It remains to see that M' n contains the uniform distribution on "simple" classes of rotation in E n , 
since any measure in A4 n is a mixture of these distributions. Take an element a = (ao, • • • , ctn-i) 
in E n . For any i S [1, nj, let a(i) = (at, . . . , a n -i, a 0i ■ ■ ■ > o>-i-i) and i? := {a(i),i G {1, . . . , n}}, 
the rotation class of a. We now show that the probability measure X^ei? $p belongs to Ai' n , 

which completes the proof. For this, denote by z% the number whose binary expansion is a(i), and 
define Q 1 and Q° to be 2 n x 2 n matr ices where all entries are except for 

It is then simple to check that fj,(x) = Trace(f3™ =1 Q Xi ) coincides with YIr^r^P- ^ 

2.2 The fundamental Lemmas 



The following two lemmas [8] and 10 - that are among the main contributions of the present 



paper - serve two different purposes. The first Lemma can be used to identify the solution //(x) = 
Trace(f3 i Q Xi ) (and/or to prove that such a measure is a solution) of (j, = f/T n when T ra has 
a product form. It gives a sufficient condition on {Q^iQ 1 ) in the form of a "finite" system of 
quadratic equations. In the generic case, this system of equations does not depend on n, and then 
a right pair (Q°, Q 1 ) will provide a representation of [i for any n. 

The second Lemma can be used when no such solution has been found. It permits us to describe 
the measure vT K on the /cth line starting from u, for some particular v using some matrices Q? re y 
Since vT K — > /i the solution of fi = fiT, this provides a way to approach \x. The problem is 

re— >+oo 

that the matrices appear to grow with k. The discussion of the convergence of the sequence 
(Q x K yK) is addressed later in this paper: this gives the clues we mentioned in the abstract that 
GF(x,y) should have an expression using some eigenvectors of some explicit infinite matrices. 
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Lemmas [8] and 10 are stated in the case of Sq(re) for processes with values in {0, 1}. They will 
adapted to the triangular lattice and to processes with more than 2 values in Section [6. 1| 

A matrix T n indexed by E n is called a transfer matrix. If for all a E E n , ^ b T ra (a, b) = 1, 
then T n is a probability transfer matrix. Moreover if, for any a, b E E n , 

n-l 

T n (a,b) = n T <^®i^ (13) 

i=0 

for some T a ^^ c , a, b, c, E {0, 1}, we say that the matrix transfer has a product form. 

Lemma 8 Assume that T n is a probability transfer matrix with a product form. Assume that there 
exists square matrices (V x , H x , x E {0, 1}) such that, V x H y = [o] for x ^ y, and such that for any 

x 

V*H*= Yl HV V y 'T y , y ,, X . (14) 

M'e{o,i} 

Finally let Q x = V X H X for x E {0, 1}, and /x(x) = Trace (n?^ 1 Q Xi ) ■ We have p = pT n . 



Note 9 (1) Equation (14) is the system of quadratic equations mentioned in the introduction and 
in the Abstract. To find a solution to fi = fiT n it suffices to find (V x ,H x ,x E {0,1}) satisfying 



(2) One may weaken the conditions in Lemma\£^ replacing (14) by "there exists an invertible matrix 
P", such that 

V x H x = p[ H y V y 'T y>yf>x J P- 1 . (15) 

V,y'£{o,i} / 

Before proving this Lemma we make two remarks. Firstly, we are only interested in non zero 
solutions! Under the hypothesis of this Lemma, the positivity of the measure fi is not guaranteed, 
nor the fact that it is non-zero or real. In the case where a non trivial solution exists, in the sense 
where /x(x) 7^ for some x E E n , then \x is a multiple of an eigenvector of T n with eigenvalue 
1. Since T n is a Markov probability kernel, by the Perron- Frobenius theorem there is exactly one 
such eigenvector, which has non-negative entries. Note that if (V x , H x , x E {0, 1}) is a solution of 
(14) then so is (cV x ,H x ,x E {0,1}) for any c / 0. Ac (which may depend on re) can be taken 
such that fj, is exactly the invariant distribution. Hence, a condition for non triviality (CNT) for a 
solution of ( 14 ) is as follows: 

CNTi : 1 is an eigenvalue of ^ V X H X . (16) 

zeo.i 



Secondly, finding a (V x ,H x ,x E {0,1}) solution of (14) is difficult in general. The number of 
variables and of equations quadratically increase with the size of the matrices (V x , H x , x E {0, 1}). 



Proof of Lemma 



Write Trace FIS, 1 = Trace nr=o Ev^e(o,H 



Since V Vi H Vi = if y\ 7^ yi, we impose y[ = Vi- It remains to write, by some usual commutations 

(n-l \ n-l 

II HViVVmi II T vumei,*i = E ^(y) T «(y' x ) = ^n- □ 
i=0 ) t=0 yeE n 

We now state the second fundamental lemma. 
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Lemma 10 Let T n be a transfer matrix on E n . Consider two matrices QJ Q \,x G {0, 1} ; of size 
m x m and some rectangular matrices VZ\,x G {0, 1} of size m x I, and Hf yX G {0, 1} of size 
I x m, such that 

Q(0) = V (0) H (0), and, for x + y, VfofZ* 0) = [0]. 

Consider the measure /i( ) defined on E n by /i( )( x ) = Trace ^rii=i Q(o)) • Assume that there exist 
8 matrices (h xy , (x, y) G {0, l} 2 ) and (v xy ,(x,y) G {0, l} 2 ) (of size i X j and j x i so that the 
products h xy v x i y i and v xy h x i y i are defined) such that for any y, x in E n , 

Trace ^11 / htt.*i<Bi%ei>*iei j = T n(y> x )- ( 17 ) 

Now, set for x G {0, 1} 

v m= E H ( )®K*, flfi)= E Qh = v m H h ^ 

ye{o,i} j/e{o,i} 
and /i(i)(x) = Trace ^ni^o 1 ■ ^ en un der these conditions, 

1) /in) = /^( )T n . If moreover fj,r \ is a probability measure and T n is a probability transfer matrix, 
then is also a probability measure. 

2) moreover if for any x,x', 

h xy v x r yl = [0] for y / y', (19) 

f/iera V/^iifL = [O] for x / y. 



5j Assume that the probability transfer matrix T ra /ias i/ie product form (13) and that there exist 
h XtV with one line, and v x>y with one column such that 

h y , x v y i tX = T y>y ' iX , for any y, y, x. (20) 



Then condition (11) is satisfied, and Q?y c,nd Qfy have the same size. If moreover Q^s = Q(q\ 



for x G {0, 1}, then /x (1) = /i (0) . 



Important note 

• The size of the cylinder n plays no role at all in the hypothesis. Hence, the same matrices 

^U) ' ^(i)' ^ or no ^ wor ^ ^ or a ^ n ' 

• In the case where = Qf Q yX G {0, 1}, we say that Q( ) is a fixed point for the matrix 

equation. That is, for any n, \x n defined on E n by 

Hn(x) = Trace y (jfaj , 

is a solution of \i n = fi n T n . 

Clearly, condition (2) is the condition needed to iterate the construction, that is to be able to 
give a matrix representation of the measure on the «'th line, as stated in the next Corollary. 
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Corollary 11 Let fJ,^)> Qfo)> ^(o)> an< ^ ^(o) sa ^ s fy the hypothesis of Lemma 10 (for example 
V/L = ff/QX := [l] and = H? ^ := [0]), and h XjV and v XiV satisfy (11), (19) and (20). Then, 
for any k > 1, x £ {0, 1}, let 



H M = ^{0,1} VI x) h 



(21) 



and Q x {k) := V^H^, x G {0, 1}, and 



'n-\ 



A t ( K )( x ) = Trace ( JJ 



(22) 



vi=0 



T/ien, for any k>\, fi^ = M(k-i)T = /^(o)T K . 



Proof of Lemma 10 , This proof relies on ( 11 ) and ( 12 ); we also make heavy use of Trace(^4i . . . Aj) 
Trace(A,-j4i . . . Aj—x). For (1), write 



Trace 

\i=l / 



(o) ^ V y'^i 



Trace (nfE^o)®^) (E< 

\»=i \ y I \ y' 

Trace LQ E V (o) ® V,** E #(Q) V*©i 

\i=l \yi J \ y 



Trace ( IJ E V (o) ® ) (#(0) ® V^i©i 

V<=1 V 2/' 



ill) ^ %>*i J (^(0) ® hvi&iG>l 



E Trace (j] (vg 
y \i=l 

/n-1 

E^(0)( y ) TraCe n 



\i=0 
'n-1 



(1) now follows from (17) 



For (2), write for x ^ y, 



\/X fry 



E^(°)(y) Trace II 'Hft«ei%ei«ei 
y \i=0 



which is indeed [0] , since h WjX v Zjy = [0] for any w, z. (3) follows immediately. □ 
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2.3 Relation with ASEP/PASEP/TASEP 



The system (14) is quite close to the systems associated with the ASEP/PASEP/TASEP prob- 
lems which use the "matrix ansatz" proposed by B. Derrida & al. [8]. In this seminal paper, it is 
shown that the invariant distribution of the ASEP can be expressed as the formal result of a compu- 
tation in some algebraic structure, where some operators D and E satisfy some quadratic relations 
of the form DE = ED + E + D (or more generally with some multiplicative parameters added) 
and some additional "border conditions" of the type D\V>= (3\V>, <W\E = a~ l <W\. For these 
problems, matrix solutions D, E, V, W are explicitly found. These solutions have, depending on the 
values of the parameters and border conditions, a finite or an infinite size. In subsequent studies 
of these exclusion processes, finite/infinite matrices appear as solutions of quadratic equations of a 
type appearing in [5]. Each time, the questions of commutation of matrices, tensor products, and 
existence of limits with some growing matrices arise. We send the interested reader to [3] and |12| . 
and references therein. 

These quadratic problems/algebras are also at the core of several different problems of enumer- 
ative combinatorics, as observed by Viennot (this is discussed at length in several of his talks and 
courses, available on his web page, see e.g. [T71 section 7]). 



3 Computation of fixed point solutions for gases of type 1 and 2 

As explained in the previous subsection, two main cases emerge. First, there may exist a fixed 



point for the matrix equation involving finite matrices according to Lemma 10 3). If no such fixed 



point solution exists (or even if it does), Corollary 11 allows one to build bigger and bigger matrices 
to describe the distribution on the kth. line. In some cases, taking the limit give rise to some infinite 
matrices. We first discuss, what happens when matrices of finite size are fixed points of the matrix 
equation. This is the case for the gas of type 1 on Sq(ra). 

3.1 Gas of type 1 

This happens in the case Z = Xs q ( n ), for which numerous different computations of the gas 
density exist. Let us add one way to find the solution. 

3.1.1 Application of Lemma [8] 



It suffices to search for finite matrices (V x ,H x ,x G {0, 1}) which solves (14). Take 

81,2 



V° 
V (0) 



S 2 ,2 



(o) - 



81,1 

s 2 ,i 



(0) - 





*1,2 *2,2 



(o) ■= 



*1,1 *1,2 





(23) 



They are chosen to trivially satisfy VZsHjU = [o] if x ^ y. 

Note 12 The choice of the letters V and H comes from these vertical and horizontal structures. 
The system of equations V X H X = V\, „/ H y V y T^,, is equivalent to: 



Sysf 



Sl,l*2,l = 0, S 2 ,ltl,l = 0, *l,iai,l = 0, S 2 ,2*2,2 = (1 - p)(si,2*l,2 + ^2,2*2,2), 
*2,lS2,l = P(S1,2*1,2 + 82,2*2,2), Sl )2 *l )2 = *1,1S1,1 + *2,lS2,l, 
Si, 2*2,2 = *l,lSl,2 + *2,lS2,2, $2,2*1,2 = *l,2Sl,l + *2,2S2,1 
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This system has non trivial solutions. For example si^ = 0,81,2 = p, S2,i = 1/(1— p), 82,2 = 1 , = 
0, ii,2 = 1) £2,1 = (1 — p)p, ^2,2 = 1 — p, in which case 



p) 




"00" 


p 




p 



Q° 



From here the density Trace((Q° + Q 1 )" -1 ^ 1 )/ Trace((Q° + Q 1 )™) of the gas of type 1 on the 
cylinder can be computed, and from that GF Sq '"'(x, 1) thanks to Proposition [3| Taking the limit 
gives GF Sq (i, 1). This representation was known and is present under different forms in [9] 13]. 
What is remarkable is that it follows from a simple computation. 



3.1.2 Application of Lemma 10 



For the same result, one may use Lemma 10 since a solution to both (20) and (19) exists. The 



idea is to search for solutions under the following form, for which ( |19[ ) automatically holds : 

h ft =[0,01,0,02], =[0,61,0,62] 
ho,i = [ci,0,c 2 ,0] , hi,i = [di,0,d 2 ,0] 



(24) 



and v xy = h xy for all x, y, this last condition not being needed. The system (20) can be rewritten 
as : 

d\ + d\ = 0, cidi + c 2 d 2 = 0, b\ + b\ - 1 = 0, 
c i + c 2 ~~ P = 0' biai + 6202 — 1 = 0, a 2 + a 2 — 1 + p = 0. 
Clearly, d\ = 0, d 2 = 0, and the rest of the system becomes 



SysJf 



(25) 



Sys x := {6? + b\ 



1 = 0, ci 2 + c 2 2 



P 



0, ai 2 + a 2 2 - 1 +P = 0, Mi + 6 2 a 2 -1 = 0}. 



Solutions to this system exist. We then take V^, Hf^ as in (23) for example, and take VZyHZx 
as defined in (21). Again, = is possible since h x ,y ana v x>y are respectively single line and 
single column (according to Lemma 10). Of course, since hy jX Vy> jX = T x (y, y', x), we are back to 



the previous problem treated in Section 3.1.1 



Note 13 The appearance of non real numbers in the considerations (Sys x has no real solution) 
does not harm the reasoning at all since T Y is still a probability transfer matrix. The possibility 
to use non real matrices V X ,H X as a solution of the equations of interest enriches the space of 
solutions. 



3.2 No finite solution for the gas of type 2 

In the case of the gas of type 2, we were not able to find finite (non-trivial) matrices (V x , H x , x G 
{0, 1}) which solve ( |14[ ). If they exist, they must have size > 5 (this can be seen via the computation 
of a Grobner basis). Also, numerical computations - which in principle does not guarantee any 
result - indicate that no non-trivial matrices (V x , H x , x G {0, 1}) with complex coefficients having 
size < 8 are solution. Hence we were unable to use Lemma [8] to go on. 



Question 1 (a) There it exist finite square matrices V°, H°, V 1 , H 1 which solve (14) forT = T Y ? 

(and such that ^2 X V X H X has eigenvalue 1?) 

(b) If not, is it possible to find a solution to (15) for T = T Y ? 
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4 Toward an infinite size solution 



This section prospective explains how some infinite matrices may arise. Even if no complete 
solution is provided, we hope that this new point of view will allow some reader to tackle the 
problem of computing GF(x,y). 

Recall that Lemma [10] and Corollary [TT] say that if there exist single line and single column 
matrices h x>y and v x>y solving h x ^ y v x ^ y r = T z (x, x', y)l y=y ', then on a cylinder of the square lattice 
the distribution on the rc+l'th line pLi R \ := fir K -i)T z has a representation of the form (22 ), provided 



that the distribution on the first line //( ) nas this same form with V^HJqs = for x ^ y. Such 
matrices h x v and v x ^ y ' exist for Z E {X, Y}: 

Theorem 14 For any Z G {X, Y}, 

(1) there exist single line matrices (h Xt y,x,y G {0,1}) and single column matrices (v Xt y,x,y E 



{0, 1}) which solves (19) and Mud, namely 



(26) 



T z (x, x', y)l y = y > for any (x, y, x' , y') £ {0, l} 4 . 

(2) For any k, the entries of the matrices Vfcy H? K y Q? K -\ ( as introduced in Corollary 11) can be 
computed. 



(3) The solutions of (26) and the initial matrices V^,H^ can be chosen in such a way that for 
any x, V(2ft)' ^(2/«)> Q(2«) converges simply when k — > +oo (that is, each fixed entry converges). 

We have already discussed the existence of solutions to Sys^, this implying the existence of 



solutions to (26) in the case Z = X. To prove Theorem [14^1), let us write the following system of 
equations for h X y and v XtV defined as in (24), in the case Z = Y. In this case (26) is equivalent to 



2 - 1, (ii 2 + d 2 2 - v + VQ - Q = 



Sys 2 



ci + c 2 — q + pq = 0, c\ + qd\ - qc\ 
d\C\ + d 2 c 2 - q + pq = 0, a\ 2 + a 2 2 - 1 + 
b\a\ + b 2 a 2 - 1 + q - pq = 0, b\ 2 + b 2 - 1 + q + p - pq = 0. 



il t- "2 -p + pq- q 
pq = 0, 



(27) 



It is not difficult to check that this system has a solution (using Maple or Mathematica, for example, 
or the computation of a Grdbner basis). 

The proofs of Theorem |14| (2) and (3) are more delicate. Their respective proofs are the object 
of Sections 14.11 and 14.21 below. 



Notation. For any pair of matrices (A , A 1 ), A* := A + A . Similarly, for any doubly indexed 
quantity a x , y , a x ^ = a Xy0 + a Xi \. 



4.1 Computations of the entries of Q^yV^lyHf k ^ 



Let Z G {X,Y} be fixed, and let h x ^ y and v XjV as defined in (24) be solution of Sys 2 (the 



sizes of h and v are 1 x m and m x 1 respectively for some m > 1). Let us compute the entries of 
^( k )i^-(k))i k — -0 Parting from some matrices V^, Hj- Q y and such that V^H^ = 
[0] for x ^ x' of size mo x tuq (for example, = HL^ = [1], V"£n = H9 Q ^ = [0]). 
For k > 2, by associativity of the Kronecker product, 



%) = E E^-2) ® v ^y) ® V* = E 



2) 



(28) 
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with 

y 

which entails, iteratively that 

V (2k)= E ^(0)^ V ^i®^^---^ V ^-i^- (29) 

ZQ,Zl,...,Z K -l 

(The matrices v z ,y have size m x m.) A similar formula exists for H? 2K y obtained by replacing h by 
v and vice versa in the previous considerations, leading to the definition of h Z)X = ^2 y h Z)V <g) v VjX . 
In the same manner, = V^H^ can be computed: write 



Using the structure of V* K _ 2 ^ and H*' K _ 2 y we get 

^(«) = E ^-2) ® <W 

where 

Qtu,x = ( v w,yh w ,z) ® {hy,x v z,x) = ^y,z,x ( v w,yh w , z 



y z y z 



(we used here that u 0j j,, ^ a> 6 is a solution of Sysf ). Again, using the same methods above 

Q (2k)= E Qlo)®^^®^,®---®^^. (30) 

Zo,Zl,...,Z K -l 

With these formulas, the entries of V^ K y ^(2k) as we ^ as °^ Q(2k) can ' 3e computed (and with 
a simple adaptation those of odd indices also). For this, recall that if 

A = B®C 

where C is a c x c matrix, B a 6 x b matrix then, for any i,j £ {0, . . . , be — 1}, 

A[i, j] = B[i div m, j div m] x C[imodm,j modm] (31) 

with the convention that for any matrix M, M[i,j] := Mj + i J+ i, and as usual xdivm and xmodra 
denote the quotient and the remainder in the division of i by m. 
Assume now that some matrices W^ 2k ^,x 6 {0, 1}, k > 0, satisfy 

W (2.)= E W^w^^w^g-^w,,.^ (32) 

Zo,Zl,...,Z K -l 

with w 1)9 having size m x m, and having size mo x mo- Therefore W^ 2/c ) has size mo x m K ; 

any z, j in {0, 1, ... , mo x m K } can be written under the following form : 

K K 

i = a K+1 (i)m K + ^m l ~ 1 ai{i), j = a K+1 (j)m K + ^m /_1 a/(j), 
i=i i=i 
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where < a K+ i(i), a K +i(j) < mo, and < ai(i),ai(j) < m for I 6 [1, k\ (apart from a K+ \ which 



may play a special role if mo 7^ m, the a/(i)'s are the digits of i in base m). Therefore, from (31), 

x w ZliZ2 [a«(i), OrO')] x • • • x w 2ft _ lA [ai(i), (33) 



ZO 2lt-l 



There is also a way to represent this with matrices, very similar to the representation of the 
distribution of a Markov chain; for any (a, b) in {0, . . . , m — l} 2 , let M w (a, 6) be the 2x2 matrix 
defined by: 

w ,o[o,6] wo,i[a,6] 
_w 1(0 [o,6] wi ) i[o,6]_ 



M w (a,6) := 



and for a, 6 £ {0, ... , mo — 1}, 



pw{a,b}:= fwj 0) [a,6] W^[a,6] 



(0)1 



We have 



W (2 K ) . J'] = Pw K+i («) , ( j)] M w [o«(i) , a«(j)] . . . M w [01 (i) , 01 (j)] 



lz=0 
lz=l 



(34) 



This fact, together with (30) and (29), proves Theorem 14 2). 



4.2 Convergence of the entries of , Hf 2l ^ , and Q( 2 k) 

We continue from the previous section. Notice that for a fixed ai(i) and all ai(j) are zero 

for large We then immediately have: 

Lemma 15 Wf 2 j] converges when k — > +00 in C /or any 2, j in N «/ and onZy if 

p w [0,0]M w [0,0]' 

converges when I goes to +00; a sufficient condition is the convergence o/M w [0, 0] . 



Again, the convergence stated in Lemma 15 will be only interesting if the limit is not zero. We 



examine the simple convergence of Vjo K y H(2k) anc ^ Q%.k) wnen K ~~ ^ +°°) to some infinite matrices 
iy^i -ff^j Q%o), meaning that, for x G {0, 1} and any i, j > 0, 

V$ K) [i,j]^VZ[i,j], Hf 2K) [i,j]^H^[i,j] and Qfafi, j] -+ QUhj] (35) 

starting with some suitable matrices ((V^, Hf Q y Q?o))' x e {0>1})- If such a convergence holds, 
the limiting infinite matrices and are moreover a solution of the following rewriting like 

1 

for xG {0,1}, Q^ = E^®<K- ( 36 ) 

w=0 



system: 



Rewriting rules such as (36) rely entirely on the corners of the matrices (Q^JO, 0], x G {0,1}). 



Writing pa/[0, 0] := [M°[0, 0], M^O, 0]], formula (36) allows us to compute the only possible corners: 

PQ^ = P Qoo M q [0,0] 
Pv x = )^Mv[0,0] 
PH^ = ^M h [0,0]. 

We examine separately the two cases Z = X and Z = Y in the next subsections. 
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4.2.1 Case of the gas of type 1. 



We work with h defined in ( 24 ) . We then find 



M v (0,0) 



c\di 
d\ 



M q (0,0) 



c\ 
d\ 



Therefore, the convergence of the sequence (VX^ k > 0) to a non zero limit is equivalent to 
d\ = 1, in which case, for any I > 1, 



(37) 



M v (0,0)' = M v (0,0) 
and the convergence of M q (0, 0) arises if c 2 = 1, in which case, for any Z > 1, 

M q (0,0)' = M q (0,0). (38) 
We may wonder if Sys^ still has some solutions if we add these conditions. 

■ The answer is yes for the condition c\ = 1 (for example, a\ = —p,d2 = 1, &i = 0,62 = l>ci = 
1, c§ = P — 1) d\ = d 2 = 0). Then it is possible to have simple convergence for Q%k} ■ 



The answer is no for the condition d\ = 1. In order to find satisfying (26) (and Z = X) 



and such that V(o,a simply converges, it suffices to increase the size of the matrices h and v defined 
in (24). Take instead 



h ,o =[0,01,0,02,0,03], hi t o =[0,61,0,62,0,63] 
^0,1 = [ci,0,C2,0,C3,0] , hi t i = [di,0,d 2 ,0,d 3 ,0] 



(39) 



and again v x>y = t h Xj y The values of M q (0,0) and M v (0, 0) are unchanged, but this time there 



are some solutions for (17) and (19) and where di = 1. Again, this is not difficult to check with 
a program like Maple or Mathematica. Note also that if we want to solve Sys^ together with the 
two equations c\ = 1 and di = 1, there exists solutions for h and v having size 6. For example: 

a\ + p = 0, a 2 = 0, a 3 = 1, 61 = 0, 6 2 = 0, 63 = 1, c\ = 1, c\ = p,c 3 = i, di = 1, d 2 = 0, d 3 = i. 

It remains to specify /? q and p v , namely the starting condition of the construction. We may 
take mo = 1 (that is, starting with lxl matrices). Taking = = H^ ' = [0] and 

QW = = = [l] leads to 



Pv 



PH = PQ ■= [0, l] • 



(40) 



With this convention, by (|37|) and (|38|), the upper-left corner of Q% K ) coincides with Q? 2re _ 2 )> an d 



the same thing hold for V and H as well. We then have pq x = pn, 



PV X 



PQ- 



4.2.2 Case of the gas of type 2. 



First, for h and v defined in (|24|), 
M q (0,0) 



{-I + p){-l + q)c\ (q + p-pq)cp 
(-1 + p){-\ + g)di 2 {q+p-pq)d\ 



M v (0,0) 



dici 
di 2 



The convergence of M q (0, 0)' to a non zero limit happens if (p— pq + q)d\ 2 + (1 — q— p+pq)^ = 1, in 
which case M q (0, 0)' = M q (0, 0) for any I > 1, and the convergence of M v (0, 0)' to a non zero limit 
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arises if d\ = 1, in which case M v (0, 0)' 



ci 
1 



for I > 1. In this case, solutions h x v of size 1x4 



exist : there exists solutions to Sys Y with the additional condition (p— pq+q)d\ 2 + {\ — q— p+pq)c 2 = 
1 or d% = 1. Again, if if we want to solve Sys Y with both conditions together there exist solutions 
for h and v having size 6, for example: 



ai 



(>2 



(l-q)(l-p),d 2 2 
p(l -q + pq) 



(l-g)(l-p), 

1 - q + pq 



(i-p)(i- q y a2 - ^- h 

i -q + pq 2 p(i-q + pq) 



do 



c 3 



(l-p)(l-g) 



6i = 0, 6 2 
l,d 3 = 0. 



0,ci = 1, 



This suffices to imply the simple convergence of Qf 2K \ > ^(Ik) ' ^(2«) " Here pQ a 
q),p -pq + q], Pv^ = [0, 1] and p Hoo = [0, 1]. 



[(-l+p)(-l + 



This ends the proof of Theorem 14 3). 



4.3 Trace of the limit and limit of the trace 



As said above, Theorem 
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gives a representation of starting from some simple (this 
could be useful to make advances on enumeration issues concerning DA with height re). The 
important and natural question is the following one : do we have, for any x = (xi, . . . , sc&), 



fn-\ 



'n-\ 



M(2«)( x ) = Trace Q*' ► ^(x) = Trace Q c 



(41) 



Since the simple convergence of Q^ 2k ) to Qoo does not imply the simple convergence of Q(2k)^(2k) 
to Q^Q^o, the answer to this question is certainly not an immediate issue. For the TASEP, the 
choices of matrices D, E satisfying the different matrix ansatz, leads or not to the convergence of 
the product (see discussions in [3(3]). 

Moreover, simple convergence of a sequence of matrices A n to some matrice does imply the 
convergence of the trace, since the trace involves an infinite number of entries. Nevertheless, P( 2k ) 
converges when re goes to +oo by Lemmajm Let /i°°(x) = lim K /i( 2 re)( x ) the limit of the measure. The 
question is: do we have p°° = p^? Since p°° is the only non trivial probability measure fixed point of 
p°° = p°°T Y , it suffices to show that p^ satisfies the same property, which would imply p^ = Xp°°, 
for some A (which must be shown to be ^ 0). In fact, by construction q yx is associated with two- row 
transitions, since Trace(n™ = i Okm) = T.^E n \i T Z,y l ®i^ T 7 l ,z ifl ,x l = (T Y ) 2 ( x ,y)- Clearly (T Y ) 2 
is also a probability transfer matrix, corresponding to an aperiodic irreducible Markov chain on a 
finite state space. Subsequently, by uniqueness, it is easy to check that p^ is the unique solution 
to n = p(T Y ) 2 . 

Using that Qoo is solution of the rewriting system ( 36 ) , if one ignores convergence and commu- 
tation issues, then 



f ii—i 



'n-\ 



f ii—l 



Trace \J[Q%\ = £ Trace \\{QI%) Trace J] q m>Xi 



.i=0 



vi=0 



.2 = 



£ Trace ( JJ Q% ) (T Y ) 2 (y,x) 

ye£„ \i=o / 



(42) 
(43) 
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and one sees that we just have to justify convergence of ni^o 1 Qoo an d the validity of the rearrange- 
ments in the infinite sum. 

We were unable to prove the validity of this. Besides, the entries of Qoo seems no to converge 
to (due to the choice of the value of M q (0, 0), needed to have the convergence of Qf K ) to Q^); 
also, seen as series in p, q, the degrees of the entries Q* ■ do not go to +00 with i,j — > 00. 

Question 2 Is it possible to prove that Trace ( ni^Q 1 Qoo ) ^ s we ^ defined (for some notion of 



convergence) and solution of (42)? 



We here review some properties of the matrices V^,H?s and of H(2k)- These properties lead 
to some questions about the structure of Q%^, and its eigenvectors (if any). First, we have 



' n—\ \ /n—1 



„ (<0 (x) = Trace UW) =fa 



v i=0 / \i=0 



This representation is very close to the standard representation of Markov chain, where here 
H?.V^ plays the role of a probability transition and 

P W : = E JT (°") V to 

a.b 

plays the role of the transition matrix. 

Proposition 16 For any k, Q*^) ^ as e W enva ^ ues 1 an d 0. The eigenvalue 1 has multiplicity 1. 

Proof. We first claim that (<^ 2k) ) 2k (1 - Q\ 2k) ) = (<2* 2k+ i)) 2k (1 - Q\ 2k+ i)) = [°] • The claim 
implies that the minimal polynomial of Q* 2k+£ ^ (with e £ {0, 1}) divides x 2k (1 — x), which implies 
that the eigenvalues are and 1. Since Trace (<3* K )) = Trace(<5* K „ 2 )), and since Trace(Q^) = 
Trace(Q* 1 ^) = 1, the eigenvalue 1 has multiplicity 1. It remains to show the claim. For this write 



\y,% J \y',x' 



) ^ V y',x> 

J \y',x' 

(k-i) (k-i) /-^ y>y ' x 

y,y' x 
and since this last sum is 1, we have 

QS0=P(^D- (44) 



Also, using ( 18 ) 



P(«) = YjQ\k-i)®Pv 



y 

where p y = v y ^h y ^. It turns out, that h yir v x ± = 1 for all x and y; hence, any product of the form 
%i,*^:ei,*%2,*^2/2,* • • • ) v y n ,*hy n ,* equals v yiik hy nik . Henceforth, k > 0, and m > 0, 

m+2 

8) Vx+flyJ. 



^ \ m+z . r ^ * m x 1 



x,x' 
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{Q\ K) ) m - 1 then {Q\ K+2) ) m+2 



Hence, if for some k > 0, and m > 0, (Q* K ))" 

The initial conditions being Q* Q ^ = [l] (and Q*-^ = [l]) we get {Q^) 1 = (Q(o))°' an< ^ then 
_ fn * ^2k fn * ^2k+i = (Q* 2k+1) ) 2k as well). □ 



'(0) 



Denote by Lr K \ and the left and right eigenvectors of Q7 K \ associated with the eigenvalue 



1. Since (Q* K )) m converges to R^L^ when 
large enough, 

(Q:r 



m 



(«) 

+oo, and since (Q* K )) m+1 



{n)> 



for m 



(45) 



Moreover Lr K \ and -R( K ) can be normalised such that Lr K \R/ K \ = 1. Notice in (45) the equality; in 
similar situations only convergence holds. For m large enough the quantity of interest 

Trace(Q[ K) (Qf K) r) = Trace(Q; K) i? {K) L (K) ). 



Using (44) 



which leads to 



pm 



(46) 



and then (since all matrices have rank 1 and Lr K \ = t Rr K )) 

L/ 



n (K-l) n (K-2) ■ 



h: 



(o)- 



Since Trace(QKi?( K) L( K) ) = Trace(#K-R( K )L (ft ) V, 1 ,) we have also 



(k)-^(k) 



H \k) H {k-1) ■ ■ ■ H \ 



(0) 



(47) 
(48) 



a triangular product whose computation seems to be quite difficult. We may also note the following 



Lemma 17 Let d (n) (l) = P(Y,, Sq(n) = 1), the density of the gas process of type 2 (this is GF Sq(n) (x, y)) 
up to change of variables, by Proposition^. For n large enough 

Ei> [1, 2i + !]!?(«) [2i + 1,1] 



Y,j>o L {K)\hj]R{K)\jA] 



RL. 



Proof. For short, we don't write the indices (k). For n large enough (by (45 )) we have (Q* 
Hence, 

d(")(l) = Trace(Qi(Q*) n_1 )/Trace((Q*) n ) = TraceiLV 1 H l R) / Trace(Lfl). 

Introduce L(l) = [L[l, mo d 2=1] and -R(l) = [R[i, l]lj mod 2=1] > the vectors L and i?, where the 
even entries are sent to 0. Now, clearly LV 1 = (LV*)(1) and ^R = (H*R)(1), then (LF*)(1) = 
L(l) and (H*R)(1) = R(l). □ 

Let us come back to the matrices and solution of the rewriting like system (36). 

Again, the value of the corner of Q x [0, 0] is given by pn[0, 0]M q [0, 0] x ~° 
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Conjecture 1 The matrix has eigenvalue 1 with multiplicity 1. Let L°° and R°° be the left 
and right eigenvector, such that L°° = t R°°. We have 

F(Y X = 1) = Y, L°°[l, 2i + 1]R°° [1, 2i + 1]/ L°° [1, j]R°° [j, 1] . (49) 

« 3 
We have seen that Qf ^ — > Qf , (simply) and Q~f has a unique eigenvalue 1, the other ones 
being 0. This convergence is not sufficient to deduce that the infinite matrix has eigenvalue 



1, and even if it is the case, the convergence of the numerator and denominator in (49) may not 
converge. In the general case, for 4 given matrices (q x , y , y) G {0, l} 2 ), the same question arises2: 
can we find the eigenvectors (and eigenvalues) of the matrix that solves the rewriting systems 



(36) 



Question 3 Is it possible to find solutions to (1^) with infinite matrices V , H , V°, H°, in the 
case T = T Y and such that moreover any product of the form Yi7=i V Xi H Xi converges ? 



5 Other similar but different considerations 



We present some alternatives to Lemma 10 Even if morally what is done has more of less the 
same taste as this Lemma, we were not able to reduce the following considerations to it. 



5.1 Research of a solution on the zigzag 

We discussed above the construction of an invariant measure relying on the research of a measure 
of the type Trace(]^[ i=1 Q x% ) on the rows of the cylinder. Two closely related constructions can be 
proposed. The first one is quite close to that discussed in [3] around the question of "Markovian 
Field". The idea is to search for an invariant measure on a zigzag on the cylinder having the 




Figure 2: Decomposition of the measure along the zig-zag. 



following forms (idea used with success by Dhar): for any (x, y) £ E^, 

n 

(m(x, y) = c n Yl d Xi>yi u yijXi@1 , (50) 
i=i 

for some complex numbers (d a> b, u a ^, a, b G {0, 1}), where "d" is chosen for "down" , and u for "up" . 
But such a measure has a cyclic Markovian structure on the lines since for m a ^ = ^ c d afi u c ^ and 
^a.b = J2 C u a,cd c ,b the induced law on the line above (up) is /x" when the measure below (down) is 
fi d with n u (x) = Cn\[m Xi ^ mi and fi d (y) = Cn\[rh yuyim . Since P(Y, = yi\Xi = Xi,X i<s i = x ie i) = 
T Xi ,x ilsl ,yi a sufficient condition for /x M = [i d and \x d = ^ ur T is that for any a, b, c, x, y 

j da,cU c ,b/ m a,b = ^a,6,o (51) 
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Notice that condition (51) can be weaken a bit since fi u = \i does not imply m xy = m xy , but 



rather rh a ^ = w a ^Tn ai b for w a ^ such that, for any x G E n (or sufficiently in the support of fi a ), 

n 

II^®i =1 - ( 52 ) 



=1 



For measures with support E n , letting iV 0) {,(x) = #{i G {0, . . .,n — 1} : (xi,Xi®i) = (a,b)}, (52) 
rewrites 

n <r (x) =i- 

(a,6)e{0,l}2 

Since on E n , N lfi = iV ,i and N 1>0 + iVi^ = #{i : Xi = 1} = n - #{i : Xi = 0} = n - N ,i - iV ,o, 
andiV^i = n — 2N\$ — Nq^, this condition is fulfilled if u>o,o = ^1,1 = 1, ^0,1^1,0 = 1> giving us 
one degree of freedom. If we deal with measures on {0, 1, 2, . . . , k}, it suffices that the product on 
all finite cycles w XljX2 . . . w Xl _ lXl equal 1 for I < k. This provides also some degrees of freedom. 



Hence, the existence of a product form as ( 50 ) is equivalent to the existence of solutions to the 
following system of equations 

m Xt yW Xt y = m X)y . (53) 
^0,0 = Wl,l = 1,^0,1^1,0 = 1 

This is a finite algebraic system, and solutions can be found using the computation of a Grobner 
basis. Again, in order to avoid trivial solutions, an additional equation has to be added: letting 
M := ( m x,y)( x , y )£{o,i}2, the equation Trace(M n ) = fi(E n ) let us sees that the existence of a non 
zero eigenvalue for M is necessary and sufficient for non triviality: 

CNT 2 : 1 is an eigenvalue of M. (54) 

The computation of the Grobner basis of the set of polynomials corresponding to the equations 



(53) gives all transitions T for which solution exists (of course, as usual T Y is not in this set). 
One can go further searching for a matrix type solution. It suffices to mix ( |51| ) with the matrix 
considerations of the previous subsection. We then search matrices (-D a ,6, U a ^, for (a, b) G {0, l} 2 ) 
of size say k x k such that for m a ^ = ^ c D afi U c ^ and rh a ,b = U a ,cD c ^-, we have 

D a ,cU c ,b = ^a,b,c m a,bj (55) 
P^^X^y — W x ,y'Wl' X ,yP 



for some k x k invertible matrix P, (w Xt y, x, y G {0, 1}) solution of (52 ). Again some non-degeneracy 
conditions must be added: define K := (mx,y)( x ,y)<={o,i} (defined by block, and having size 2k x 2k). 
The needed condition is Trace(M n ) ^ 0, for n large. By triangulation of M, it appears clearly 
that a necessary and sufficient condition is that M has a non zero eigenvalue, which amounts to 
imposing CNT2. 

We were unable to show existence/non-existence of such matrices (P, D, U) and weights (w x ^ y ) 



solution of (55) + CNT2 in the case T = T Y for m > 3. 
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5.2 Research of a solution by projection 

Recall a simple fact : if (Xi,i > 0) is a Markov chain with state spaces S (with #S > 2), and 
if (p '■ S — > S' , then in general (4>(Zi),i > 0) is not a Markov chain. Here, one may search if the 
process Y (on E n ), gas of type 2, whose measure \i is solution of /i Y = /i Y T Y , can be written 
Yi = 4>{Zi) for some function eft, and some process (Zi, 1 < i < n) taking its values in a state space 
S with \S\ > 2, and having a simple representation. We were not able to find such a solution. 



6 Extensions 

In this section we discuss some extensions of the method we developed above: first to the 
triangular lattice, second to processes with more than 2 states. 



6.1 Triangular lattice 



We proceed as in Section 5.1 where a zigzag is considered. Since no new result are provided for 
DA on the triangular lattice, we just explain how the previous method can be adapted here. First 
the two probabilistic local transitions needed for the definition of the gases of type 1 and 2 on the 
triangular lattice are 



T&. d = F(B p (l-a)(l-b)(l-c) = d) 
Tj bc4 = F(B p abc+(l-B p )B q = d). 

Again Proposition [3] says that the density of the corresponding gas provides up to a change of 
variables GF(x, 1) and GF(x,y), the area, and area-perimeter GF of DA on this lattice. As 




Figure 3: Representation of the transition. DA on the triangular lattice; the white circles represent 
perimeter sites. On the second picture, a DA on the triangular lattice with periodic conditions; on the 
last picture, a zig-zag on the triangular lattice 



explained in Section 5.1, one searches for a representation of the zigzag process distribution as 
follows: 

A*ui,di ) « a ,...,«n I *, = Trace (d u ^ . . . U d - U A (56) 



for some matrices D u,d and JJ d,u for u,d G {0, l} 2 (see Figure [3] to see the respective positions of 
the Uj's and di's). To find some finite matrices doing the job, it it sufficient to find matrices D u,d 
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and U d,u for u,d G {0, l} 2 solving the following system: 



D abjjbc = Q; jjab D bc = Q) for G) b> c g | Qj 1 ^ 

£>4ftf/tf = ^ u U du D ud 'T dud ,j for d, /, d! G {0, 1}. 



(57) 
for 



any x, y, D 



i,i 



, D ' 1 = 


—rp 


rp 


, D°'° = 


1 r 




—rp 


rp 




I r 



, where r 



For T = T x there is a solution with "matrices" D x ' y , U x,y of size 2, for example D x,y 

— 1 —r 
-(rp + 1 + 2p) jp 1 

satisfies p + (1 — 2p)r + r 2 p = 0. This leads to the searched density (similar approach are present 

in pa g). 

Again no such chance arises for T Y . To adapt the construction of growing matrices as explained 
in Corollary 11 some single line and single column matrices h a ^ c , v a & c indexed by (a, 6, c) G {0, l} 3 



which solves the following system must be found: 

The idea then is to grow the matrices Dt K \, Ur K \ as follows: 

= Ex 



(58) 




2^x> u ( 



'(«) 



hdxz 

5 



Under this condition if /i( K ) has a representation as that given in (56) with some matrices Ur K \ 
(instead of D, U) such that DfKUff = UfKD^f = [o] if y ^ y' , then this property is inherited for 
^(re+i)) ^(k+i)- O n the zigzag below (as drawn on Figure [3]), then measure /^( K+ i) will also be given 
as in (56), with Dr K+1 \, f7( K +i) instead of Z?( K ), Ur K y Again, solutions to the system (58) exist: this 



permits one to grow some matrices D^,U^, and then to have a representation of the measure 
Infinite matrices and appear again by some passage to the limit. We were not 

able to deduce from them GF. 



6.2 Bond percolation 

Let A be a DA on G = (V, E). A bond in A is an edge e£fi between elements of A; let N(A) 
be the number of bonds in A. The GF S^(x,y) = ^2 A x^ A ^y^ N ^ of DA with source C counted 
according to the area and number of bonds can be obtained also by computing the density of some 
gas process (this is explained in [JJ, page 21, in the case G = Sq(n)). Indeed, for any G G Q, 

B% = ^ £ B%y C ^ (59) 

DcCh(C) 

where C — > D is the number of bonds between C and D. For any cell d G Ch(C), let Sc{d) = 
{(c, d) G -E |c G C} be the set of bonds from C to c£. Further, denote by Dc{i) = {d £ D \ \Sc(d)\ = 
i} be the subset of D of cells being extremities of i edges coming from C. We have 

b§ = x w ^ B G y i:i\Dc(i)\_ m 

DcCh(C) 

We will define a gas whose density will coincide up to a change of variables to Be- For this, 
associate with the set of vertices V i.i.d. Bernoulli(p) random variables (denoted (B®,x G V)), and 
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with the edges of E, i.i.d. Bernoulli(g) random variables (denoted (B®,x £ E). Consider now the 
gas defined by 

x x =b; n i-Ud n B i x,d) )- ( 6i ) 

deCh(x) \ (x,d)£E J 

Taking the expectation in the previous line, leads to 

f(x x = i,xgC) = pl c l Pft = Uef)(-i) |fl| g Ei|Dc(,)| . 

DcCh(C) 

Set now Hc(p,q) = (—1)^F(X X = l,x £ C). The last equation rewrites 

H C (p,q) = (-p) lCl £ H d (p,Q)q' £iIDc(i)1 

DcCh(C) 

in other words, the family (Hc(—p,q),C) satisfies the same equations as the family (Bc(p,q),C), 
the initial conditions being B = 1 and H = 1. On the square lattice, set 



T* y>z = P(X = z\X x = x,X 2 = y)= F(z = B p (l - xB^)(l - yB, 



Again, there exist solutions to the equations h xw v x ^ y i = l y=y iT^ x , , and we may also impose that 
M q (0,0)' converges. This permits again to use Corollary [IT] to represent ^ K y 

6.3 Bicolouration 

We present a new model of gas X taking three values 1, 2 or 3 (we see the value as a colour), 
having an interest from a combinatorial point of view, and illustrating the universality of the present 
approach. Let G = (V, E) be in G, and let (C x ,x £ V) be a i.i.d. random colouring of the vertices 
of G, such that P(C = i) = p% for i < 2, and YlPi = 1- We- set 

X X :=C X J] l Xc ^ Cx . (62) 
cech(x) 

The arguments given in Lemma [2] implies that X is a.s. well defined for p\ + p 2 small enough. 
According to (62), if C x = then X x = (with proba. 1), and if C x = % for i £ {1,2}, then X x = i 
with proba. pi if no child of x has colour i, and X x = in the other cases. 

Clearly P(Xi = 1) = -GF G (-pi,l) and F(X = 2) = -GF G (-p 2 ,l) since the gases Y = 
lx=\X and Y' = lx=2X/2 have the same transitions as the gas of type 1 with parameters p\ and 
P2 (and by Lemma p5J). This gas is related to the counting of bi-coloured DA, having no bicoulor 
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neighbouring sites (this model is interesting only when several sources are involved). 

Formally consider two sets S\ := {ci, . . . , Cfc} and S2 ■= {di, . . . , di} such that SiCi S2 = and 
such that Si U S2 is a free set. We call bicoloured DA a pair (A, I), where I : A — > {1, 2} (where 1(a) 
is seen as the colour of a). A bicoloured DA is said to be well coloured, if for any (a, b) G A 2 n -E, 
/(a) = 1(b), meaning that neighbouring sites have the same colour. 

Denote by As 1: s 2 the set of well coloured DA (A, I) with source Si U S2 such that l(Si) = {i} 
(the cells of the sources Si have colour i). This model is hard to deal with using heap of pieces 
arguments. Let GF s 1: s 2 De the corresponding GF, counted well coloured DA according to their 
number of cells of each colour. 

Proposition 18 We have 

GF Sl>Sa (-pi,-p2) = (-l)* Sl+ * S2 F(X x = l,x€S u X x = 2,x€ S 2 ). (63) 

Here again, the gas transition T(a, b, c) = F(C X = c \C Cl = a, C C2 = b) can be written T(a, b, c) = 
h a ,cVb,c'^c=c' for some monoline and monocolumn h and v, for any a, b, c, d G {0, 1, 2}. Again, this 
ensures the existence of a representation of the measure on the Kth line of the lattice using some 
matrices V* K yH* K yQ* K y for x G {0, 1, 2}, starting for some measure /X( )(x) = Trace(n™=o ^(oj"^(o)) 
and some matrices V"^, -ff^ )' sucn that V^if^ = [0] for x 7^ y. 

Here the case is particularly interesting: P(Xj; = j) for j G {0, 1,2} as well as P(X X = j, x G C) 
are easy to compute: the reason is that the projection Y and Y' (as defined above) are well known, 
and have a simple product form: they correspond in the first case to identify the states 2 and 
and in the second one to identify and 1. Hence, both Y and Y' are Markovian on a line of the 
lattice (and hard particle model on the zigzag), but X is not Markovian on the lines (or on the 
zigzag). The combinatorial issue is not to find the density of X, but rather to compute a quantity 
as P(Ao = 1, Xi = 2). At this moment, I am not able to do this. 

Thanks : Grateful thanks to David Renault for numerous stimulating discussions about this work. 
Many thanks are due to the anonymous referees for their comments. 
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